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Abstract 

We present a precision study of the ip' — > ir^J/ip and r]J/ip decay modes. The measurements are 
obtained using 106 x 10 6 ip' events accumulated with the BESIII detector at the BEPCII e + e~ collider 
operating at a center-of-mass energy corresponding to the ip' mass. We obtain B(ip' — > Tr°J/ip) = 
(1.26 ± 0.02 (stat.) ± 0.03 (syst.)) x 10~ 3 and B(ip' -»• r)J/ip) = (33.75 ± 0.17 (stat.) ± 0.86 (syst.)) x 
10" 3 . The branching fraction ratio R = B ^7^°jf$ is determined to be (3.74 ± 0.06 (stat.) ± 
0.04 (syst.)) x 10 -2 . The precision of these measurements of B(ip' — > 7r°J/tp) and R represent a 
significant improvement over previously published values. 



PACS numbers: 13.25.Gv, 13.20. Gd 



I. INTRODUCTION 



The study of the hadronic transitions between charmonium states has been an active field 
both for experimental and theoretical research. The decays if)' — > rjj/ip and 7T°J/if) were first 
observed thirty years ago, and improved measurements of the corresponding branching fractions 
were performed by the BESII [l] and CLEO Q collaborations. These decays are important 
probes of if)' decay mechanisms that are characterized by the emission of a soft hadron. The 
QCD multipole-expansion (QCDME) technique was developed for applications to these heavy 
quarkonium system processes. For this, the measured branching fraction for if)' — > rjj/if) can 
be used to predict the i] transition rate between T states 0]. 

The branching- fraction ratio, R = gfez^j^y > with B denoting the individual branching 
fraction, was suggested as a reliable way to measure the light-quark mass ratio m u /md |4|. 
Based on QCDME and the axial anomaly, the ratio is calculated to be R = 0.016 with the 
conventionally accepted values of the quark masses m s = 150 MeV/c 2 , rrtd = 7.5 MeV/c 2 and 
m u = 4.2 MeV/c 2 [5j. Previously published measurements of this ratio give a significantly larger 
value of R = 0.040 ± 0.004 jfj)]. Recently, using chiral-perturbation theory, the Jiilich group 
investigated the source of charmed-meson loops in these decays as a possible explanation for 
this discrepancy Under the assumption that the charmed- meson loop mechanism saturates 
the ip' — > it°{ri)J /ip decay widths, they obtained a value R = 0.11 ± 0.06, which indicates 
that the charmed-meson loop mechanism can play an important role in explaining the data. 
With parameters introduced into the charmed-meson loop fixed using B(ip' — > T]J/if>) as input, 
the hadron-loop contribution to the isospin violation decay if)' — > 7i°J/ip can be evaluated jl, 
0]. Measurements of these branching fractions can provide experimental evidence for hadron- 
loop contributions in charmonim decays, and impose more stringent constraints on charmed- 
meson loop contributions. It will also help clarify the influence of long-distance effects in other 
charmonium decays, e.g. -?/>(3770) — > ir°(ri)J/if) 0, 13], if)' —> ^ri C) and J/if) — > 7% 11 . 

This paper presents the most precise measurement of the ratio R and the related branching 
fractions for if)' — > 7t°J/if> and r)J/if>. 



II. BESIII EXPERIMENT AND DATA SET 



The BESIII experiment at the BEPCII 12| electron-positron collider is an upgrade of BE- 



SII/BEPC [13(1 . The BESIII detector is designed to study hadron spectroscopy and T-charm 
physics [14] . The cylindrical BESIII spectrometer is composed of a Helium gas-based drift cham- 
ber (MDC), a Time-Of-Flight (TOF) system, a CsI(Tl) Electromagnetic Calorimeter (EMC) 
and a RPC-based muon identifier with a super-conducting magnet that provides a 1.0 T mag- 
netic field. The nominal detector acceptance is 93% of 47r. The expected charged-particle 
momentum resolution and photon energy resolution are 0.5% and 2.5% at 1 GeV, respectively. 
The photon energy resolution at BESIII is much better than that of BESII and comparable to 
that achieved by CLEO [l5[ and the Crystal Ball [l6|. An accurate measurement of photon 
energies enables the BESIII experiment to study physics involving photons, 7r° and rj mesons 
with high precision. 

We use a data sample of (106.41±0.86) x 10 6 if)' decays [l7| . corresponding to an integrated 
luminosity of 156.4 pb" 1 . In addition, a 43 pb _1 data sample collected at 3.65 GeV is used for 
QED background studies. 

To optimize the event selection criteria and to estimate the background, a GEANT4-based 



simulation 18] is used that includes the geometries and material of the BESIII detector com- 



ponents. An inclusive ip' decay Monte Carlo (MC) sample is generated to study backgrounds. 



The generation of if)' resonance production is simulated with the MC event generator kkmc 



20 



for known deca y m odes with branch- 



[191 ]. while if)' decays are generated with besevtgen 

ing fractions set to the world average values [fl, and with lundcharm [2jJ for the remaining 
unknown decays. The analysis is performed in the framework of the BESIII offline software 
system [22j which handles the detector calibration, event reconstruction and data storage. 



III. EVENT SELECTION 

Selection criteria described below are similar to those used in previous BES analyses 



24j. Candidate tc° and r\ mesons are reconstructed using two photons 77, and the J /if} is 
reconstructed from lepton pairs l + l~{l = e or ji). 

Photon candidates are reconstructed by clustering EMC crystal energies. The energy de- 
posited in nearby TOF counters is included to improve the reconstruction efficiency and the 
energy resolution. Showers identified as photon candidates must satisfy fiducial and shower- 
quality requirements. A minimum energy of 25 MeV is required for barrel showers (| cos 9\ < 0.8) 
and 50 MeV for endcap showers (0.86 < | cos#| < 0.92). Showers in the angular range between 
the barrel and endcap are poorly reconstructed and excluded from the analysis. To exclude 
showers generated by charged particles, a photon is required to be separated by at least 10° from 
the nearest charged track. EMC-cluster timing requirements are used to suppress electronic 
noise and energy deposits unrelated to the event. The number of photons, iV 7 , is required to 
be V 7 > 2. 

Charged tracks are reconstructed from hit patterns in the MDC. The number of charged 
tracks is required to be two with zero net charge. For each track, the polar angle 9 must satisfy 
I cos ^| < 0.93, and the track is required to originate from within ±10 cm of the interaction 
point in the beam direction and within ±1 cm of the beam line in the plane perpendicular to the 
beam. The J /if) — > l + l~ candidates are reconstructed from pairs of oppositely charged tracks. 
Tracks are identified as muons (electrons) if their E/p ratios satisfy 0.08 c < E/p < 0.22 c 
(E/p > 0.8 c), where E and p are the deposited energy in the EMC and the momentum of the 
charged track, respectively. 

To reduce the combinatorial background from uncorrelated 77 combinations and to improve 
the mass resolution, a four-constraint kinematic fit (4C-fit) is applied with the hypothesis 
if}' — > 77Z + /~ constrained to the sum of the initial e + e~ beam four-momentum. For events with 
more than two photon candidates, the combination with the smallest y 2 is retained. 

The invariant-mass distribution for lepton pairs (Mu) is shown in Fig.IU, where the J/if> signal 
is clearly seen with a high signal to background ratio. For the further analysis, events are kept 
for which the reconstructed J/if> mass falls within a window of Mu £ (3.05, 3.15) GeV/c 2 ; a 
mass window that is significantly larger than the mass resolution of about 8 MeV/c 2 . Figure [2] 
shows a Dalitz plot of the invariant-mass squared M^ hJ ^ for the reconstructed J/if> and the 
energetic photon versus the two-photon invariant-mass squared where 7^ denotes the 

photon with the higher energy E lh > E T Bands of 7r°, rj, and Xcj (J = 0,1,2) are clearly 
visible. To suppress the dominant source of background, which is from Xcj decays, the mass 
of the 'jhJ/ip system is required to satisfy the condition M lh j/^ ^ (3.50, 3.57) GeV/c 2 and 
M Jh j/ip < 3.5 GeV/c 2 for if}' — > 7r°J/if> and r]J/ip, respectively. The 4C-fit \ 2 is required to be 
less than (100, 60, 70, 50) for the final states (ii e + e~, 7r°/i + /i _ , r]e + e~, i]fi + fi~), respectively, 
where the values are determined by optimizing the statistical significance S/ a/S 1 + B, with S(B) 
the number of signal (background) events. The background event levels are determined from 
the if}' inclusive MC sample. 
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Fig. 1: The invariant-mass distributions for (a) electron-positron and (b) di-muon pairs in the selected 
77/+/" events in the data. 




Fig. 2: Dalitz plot of M 7fc7 (vertical) versus M^ h j/f (horizontal) for data, where 7^ denotes the ener- 
getic photon. The horizontal bands around M 2 h7 =0.02 (0.30) (GeV/c 2 ) 2 are due to ip' — > ir°(r]) J/ip 
transitions. The vertical bands around M^j,^=11.65 (12.30, 12.70) (GeV/c 2 ) 2 are due to transitions 
■0' — > 7XcO(ci,c2); the arrows denote the requirements to remove backgrounds from the Xc\,c2 states as 
described in the text. 

IV. DATA ANALYSIS 

Background events from ip' decays are studied using the inclusive MC sample. The back- 
ground is dominated by ip' — > ^Xcj, Xcj — > iJ/ip ~^ decays. In addition, there are a 
few events from direct ip' — > ^J/ip, J/ip — >■ decays j25|. The shape of the M 77 distri- 
bution of direct 77 J/ijj decays is smooth within both the tt° and the r] mass regions. The 
non-resonant background from ip' — > 77/+/" is studied using J/0-mass sidebands in the data. 
For ip 1 — > 7] J/0, there is an additional background from 0' — > 7r°7r°J/0, which has a smooth 
shape within the ?7-mass region. The background contribution from QED processes is studied 
using the continuum data taken at y/s = 3.65 GeV, and it is found to be negligible. The sum 
of all the MC-determined backgrounds in the M 77 distribution are shown in Fig. [3J and they 
are found to be in reasonable agreement with those observed in the data. 

To determine the detection efficiency, the angular distributions are properly modeled in the 




event generator which accounts for polarization in the ip' and J/ip decays. These decays are 
dominated by their transverse polarization; longitudinal polarization of the ip' is negligible due 
to its production from e + e~ annihilation, and, since the J /if) is produced via ip' — > n°(n)J/ip 
transitions, its longitudinal polarization vanishes because of parity conservation. Thus, their 
polar-angle distributions take the form of dN/d cos 9 oc (1 + cos 2 9), where 9 is the polar angle of 
J/ip in the ip' rest frame for ip' — > n(n°)J/ip decays, or the angle between the lepton momentum 
in J/ijj rest frame and the J/ip momentum in the ip' rest frame for J/ip — > decays. As an 
example, Fig. H] shows angular distributions for J/ip and fi~ in ip' — » r/J/ip — > r]u + u~ decays, 
where the angular distributions obtained from MC simulations (histograms) are observed to be 
in excellent agreement with the data (the points with error bars). Similarly, we have verified 
that the angular distributions in the ip' — > Tr°J/ip decay are well described by MC simulations. 
The detection efficiencies are determined using these MC event samples, and the values are 
listed in Table [B The efficiencies for / -f-fe + e~ final states are lower than that for / -f-fu + u~ final 
states because the e + /e~ tracks suffer from stronger bremsstrahlung effects. 

The signal yields are obtained from fits to the observed two-photon invariant mass M 77 
distributions. The observed line shapes are described with modified ir°/n line shapes plus 
backgrounds. The 7r° and n line shapes are taken from MC simulation including detector reso- 
lution; the 7T° and n are described with relativistic Breit- Wieners in the event generation, and 
their masses and widths are fixed at their nominal values |6j. To account for possible differ- 
ences in mass resolution between data and MC simulation, the ir°/r] line shapes [LS(7r /^)] are 
modified by convolution with a Gaussian function G(M 77 — 5m, a). This technique of mass 
resolution smearing treatment was used in a previous publication [24|. The probability distri- 
bution function (PDF) for the signal is taken as LS (7^/77) £g>G(M 77 — 5m, a), where 5m and 
a correct the n°/n mass and mass resolution, respectively. The PDF for the dominant back- 
ground contribution is obtained from MC simulation and the residual background contribution 
is modeled as a first-order and a second-order polynomial function for the 7r° and n channels, 
respectively. The polynomial coefficients are free parameters with values determined from the 
data. The fit results are shown in Fig. [31 For ip' — > 7r°J/ifj, the fit yields 1823±49 events 
for the J /ip — > e + e~ sample with a goodness of fit of x 2 /ndf = 0.85, and 2268±55 events for 
J/ip — > fi + fi~ with a x 2 /ndf = 0.86, where ndf denotes the number of degrees of freedom in 
the fit. For ip' r/J/ip, the fit yields 29598±202 events for J/ip e+e~ with a y^/ndf = 1.33 
and 38572±280 events for J/ip — > u + fi~ with a x 2 /ndf = 0.96. The resulting values of 5m and 
a are \5m\ < 1 MeV/c 2 and a < 3 MeV/c 2 for all the modes. The signal yields are listed in 
TableE 

The branching fractions are calculated from the expression 

where X represents 7r° or n, N sig and N^i are the signal yields and the number of ip' events, 
N^i = 106.41 x 10 6 . B(X — > 77) and B(J/ip — > denote the branching fractions of 

n° /n — > 77 and J/ip — »■ e + e~(/i + /i~) [fjj]. The variable e represents the MC-determined detection 
efficiency. The measured branching fractions for each final state are listed in Table [TTl 

To validate the event selection criteria and fitting procedure, we perform a study using a MC 
sample of 106 xlO 6 inclusive ip' events, with the known branching fractions as input. The same 
analysis procedure as used for the real data is applied for this MC sample and the obtained 
branching fractions for the ip' — > ir°(n)J/ip channels are found to be consistent with the input 
branching values within the statistical accuracy. 




Fig. 3: (Color online) M 77 distributions and fit results, (a) ip' — > ir°J/ip, J/ip — > e + e~, (b) ip' — > 
7r°J/ip, J/ip — > fJ. + fJ-~, (c) ip' — > rjJ/ip,J/ip — > e + e~, (d) ip' — > rjJ/ip, J/ip — > fi + , where the points 
with error bars are data, and the solid (red) curves are the total fit results, and the dashed curves 
are the fitted background shapes. The hatched histograms represents dominant background events 
obtained from MC simulation and J/ip mass sidebands. 



A. SYSTEMATIC ERRORS 



The main sources of systematic uncertainty originate from the number of ip' events, the 
trigger efficiency, the lepton tracking, photon reconstruction, kinematic fitting, uncertainties of 
the branching fractions for tt°(i]) — > 77 and J/ip — > e + e~(yU + /i~), and the selection and fitting 
procedures. 



The uncertainty on the number of ip' events is 0.81% as reported in Ref. [17|. Trigger 
efficiency uncertainty is 0.15% as reported in Ref. 26[. The photon reconstruction uncertainty 
is determined to be 1% per photon in Ref. {23}, and, thus, the two-photon final state is assigned 
an uncertainty of 2%. The tracking efficiency of the hard leptons is studied using a control 
sample of ip' — > n + n^J/ip, J/ip — )■ e + e~(/i + /i _ ) decays. The tracking efficiency e is calculated 
as e = N f u u/ N a u, where Nf u u indicates the number of tt + tt~1 + 1~ events with all final tracks 
reconstructed successfully; and N a u indicates the number of events with one or both charged 
lepton tracks successfully reconstructed in addition to the pion-pair. The difference in tracking 
efficiency between data and MC is calculated bin-by-bin over the distribution of transverse 
momentum versus the polar angle of the lepton tracks. By this method, tracking uncertainties 
are determined to be 0.14% (0.20%) and 0.16% (0.19%) for ip' ->- n°J/ip, J/ip ->- e+e-(//+/i~) 
and ip' —> T]J/ip, J/ip —> e + e~ respectively. 

Some differences are observed between data and MC x 2 distributions from the kinematic 
fit. These differences are mainly due to inconsistencies in the lepton track parameters between 
MC and data. We apply correction factors for various ^ (e ) track parameters that are ob- 
tained from control ip' — > 7T + 7T~J/ip data samples, where J/ip — > e + e~(fi + fi~). The correction 
factors are found by smearing the MC simulation output so that the pull distributions prop- 




cos9(JA|/) cos6((i ,JA|/) 

Fig. 4: (Color online) Angular distributions for (a) J/if> in the ifi' rest frame, (b) [i~ in the J /if) 
helicity system, where 0(fJ,~, J/ip) is the angle between the \x momentum in J /if) rest frame and the 
J/ip momentum in if>' rest frame. Points with error bars are data, and histograms are MC simulations 
as described in text. 

erly describe those of the experimental data. Half of the differences between the detection 
efficiencies, obtained using MC simulations with and without applying these correction factors, 
are taken as systematic errors. These are 0.15% (0.19%) and 0.20% (0.28%) for if)' — > ir J/if), 
J /if) —> e + e~ n~) and if)' —> r]J/if>, J/if> — > e + e~ respectively. 

Requirements on the E/p ratio and the invariant mass Mi+i- have been applied in the event 
selection. Uncertainties associated with these requirements are determined using the same 
control sample described above. Differences in the detection efficiency between the control data 
sample and the MC due to the E/p ratio requirement are 0.06% and 0.05% for J /if) — > e + e~ 
and J /if) — > /i + yU~, respectively. Uncertainties caused by the mass window selection are 0.06% 
for both the e + e~ and /i + yU~ channels. 

An uncertainty due to the M 7 ;+ z - requirement arises from a difference, S(x c i,2), m the x&,2 
mass resolution between the data and MC simulation, and is estimated by changing the MC- 
optimized requirement to one optimized using the data. Uncertainties caused by this M y i+i- 
requirement are determined in this way to be 0.55% (0.16%) and 0.11% (0.82%) for if)' — > TT°J/ip, 
J/ip —> e + e~ (n + n~) and ip' —> rjj/ip, J/ip — > e + e _ ( / u + / u _ ), respectively. 

Systematic errors due to the background shape are estimated by varying the function used 
to describe nondominant backgrounds from a 1st- (2nd)- order polynomial to a 2nd- (3rd-) 
order polynomial for if)' — > it (77) J / 'if) . The difference in the signal yields observed is taken as a 
systematic error. The uncertainty due to the choice of fitting range is estimated by repeating 
the fits using a fitting range that is 80% as wide as that used in the original fit. The difference 
in the signal yields is taken as a systematic error. Table II III summarizes all the sources of 
systematic uncertainties. 

B. RESULTS AND DISCUSSION 

Branching fractions for the decays if)' — > 7i°J/ip and r)J/ip with J/ip — > e + e~, fi + fi~ 
are calculated with Eq. (TjQ) using the fitting results and the detection efficiencies as inputs. 
Branching fractions measured using J /if) — > e + e~ and /i + yU~ final states are combined to- 
gether with the weighted average method described in Ref. [gj], here common systematic un- 



certainties are counted only once. The combined branching fractions are B(ip' — > n^J/ip) = 
(1.26±0.02±0.03) x KT 3 and B{^ ->• tjJ/i/)) = (33.75±0.17±0.86) x 10~ 3 (see TableE]). Using 
the measured branching fractions, the ratio R = B(ip' —> n°Jip)/B(tp' — > riJip) is calculated to 
be R = (3.74 ± 0.06 ± 0.04) x 10~ 2 (see Table HT1) . Note that systematic uncertainties that are 
common to both channels cancel in the ratio. Our combined result on the i?-ratio is consistent 
with previous world average values with a precision improvement of about a factor of five. 

These precise measurements of the ip' — > 7r°J/?/> and rjJ/ip branching fractions permit the 
study of isospin violation mechanisms in the ip' — > -K Q J/ip transition. As shown in @, the ax- 
ial anomaly does not adequately explain the observed isospin violation, while contributions from 
charmed meson loops would be a possible mechanism for additional isospin violation sources. 
Confirmation of sizeable contributions from charmed-meson loops would be an indication that 
non-perturbative effects play an important role in the charmonium energy region. 

Table I: Summary of signal yields and detection efficiencies for the each final state. 



Mode 


lj/ -»■ TT° J/lp 




rjj '/ij) 


Final state 


77e + e~ 77/i + /i~ 


77e + e~ 


77,u + M~ 


e{%) 


23.05 29.11 


35.41 


46.28 


N si 9 


1823±49 2268±55 


29598±202 


38572±280 



Table II: Summary of measured branching fractions (B) and the ratio R = ^g^/z^j^y with comparison 
to world average values (see PDG). 



B or R 


Final state 


This work 


Combined 


PDG [6] 


B{$' -»■ n° J/ip) 
(xlO" 3 ) 


77e + e~ 


1.27 ±0.03 ±0.03 
1.25 ±0.03 ±0.03 


1.26 ±0.02 ±0.03 


1.30 ±0.10 


B{ip' TfJ/ifi) 
(xlO^ 3 ) 


77e + e~ 
77/i + //~ 


33.77 ±0.23 ±0.93 
33.73 ±0.24 ±0.90 


33.75 ±0.17 ±0.86 


32.8 ±0.7 


D _ B(#->7T U J/V) 

(xlO" 2 ) 


77e + e~ 
77/i + /i~ 


3.76 ± 0.09 ± 0.06 
3.71 ±0.09 ±0.05 


3.74 ±0.06 ±0.04 


3.96 ±0.42 



Table III: Summary of all systematic errors (%) considered in this analysis. 



Sources 


7r°J/V>(e + e~ 


) ir J/i>(ti+n- 


) r]J/ip(e + e 


) 7?J/V'(/U + n ) 




0.81 


0.81 


0.81 


0.81 


Trigger 


0.15 


0.15 


0.15 


0.15 


Tracking 


0.14 


0.20 


0.16 


0.19 


F* Vint, on 


2.00 


2.00 


2.00 


2.00 


4-C Fit 


0.15 


0.19 


0.20 


0.28 




1.01 


1.01 


1.01 


1.01 


B r (ir°/ri -»• 77) 


0.03 


0.03 


0.51 


0.51 


M(i+rj 


0.06 


0.06 


0.06 


0.06 


M( 7 /+r) 


0.55 


0.16 


0.11 


0.82 


E/p 


0.06 


0.05 


0.06 


0.05 


Background shape 


0.24 


0.24 


1.14 


0.10 


Fitting range 


0.63 


0.80 


0.55 


0.58 


Total 


2.55 


2.55 


2.77 


2.66 
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